Thermodynamic Consistency of the g-Deformed Fermi-Dirac Distribution in 

Nonextensive Thermostatics 
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The 5-deformed statistics for fermions arising within the non-extensive thermostatistical formalism 
has been applied to the study of various quantum many-body systems recently. The aim of the 
present note is to point out some subtle difficulties presented by this approach in connection with the 
problem of thermodynamic consistency. Different possible ways to apply the q-deformed quantum 
distributions in a thermodynamically consistent way are considered. 



The g-deformed quantum distributions [ll-[lol| inspired 
on the nonextensive thermostatistical formalism |ll| - [l3j 
have been the focus of considerable attention in recent 
years (see p^ - [l9j and references therein). An interest- 
ing recent development along these lines was the for- 
mulation by Pereira, Silva, and Alcaniz (PSA) of a q- 
deformed equation of state for relativistic nuclear matter 
within Walecka's phenomenological relativistic approach 
[Tg} . The PSA equation of state may be relevant for the 
study of nuclear matter in neutron stars. Strictly speak- 
ing, however, the PSA equation of state as derived in [l^ 
is not thermodynamically consistent. Here we will con- 
sider possible solutions to this difficulty. As we are going 
to demonstrate, the PSA approach can be implemented 
in a thermodynamically consistent way either by adopt- 
ing a g-nonlinear form for thermodynamical quantities 
like the total energy E or the total number of particles N 
that in the standard thermostatistical formalism are lin- 
ear functions of the (mean) occupation numbers or, alter- 
natively, by recourse to a different choice of the entropy 
functional. This last procedure can, in turn, be imple- 
mented in two different ways. One possible re-definition 
of the entropy, accounting for the relevant equations of 
constraint, yields the same quantum distribution func- 
tions as solutions of the entropic variational principle as 
the one employed by PSA, and it preserves thermody- 
namic consistency. 

Consider the g-deformed Fermi-Dirac distribution used 
by PSA which has the form 



l + [l + (g-l)(a + /3e,)]^ 
1 

l + [l + (g-l)/?(e, -A*)]^' 



(1) 



where are the single particle energies, (5 = 1/kT (T 
being the absolute temperature and k denoting Boltz- 



mann's constant), /i 



/3 



is the chemical potential, and 



if a + /3 > 



(2) 



As shown in [20|, the Fermi mean occupation numbers 
(H)) can be obtained from a maximum entropy principle 
based on the entropic measure 



i 

where the function Cq{x) is defined by 



(3) 



if x< 



if X > i 



(4) 

The optimization of the entropic measure ([3]) under ap- 
propriate constraints corresponding to the total number 
of particles N and the total energy E leads to the varia- 
tional problem 



k 



N- 



\+(3\ E~ 



(5) 

where a and /3 denote the Lagrange multipliers asso- 
ciated, respectively, with the aforementioned two con- 
straints. The solution of equation ^ is given by the q- 
deformed Fermi distribution ([T]). In the limit case q 1, 
the entropic functional ([3]) reduces to the well known 
Fermi functional 



; Inrii -I- (1 - Hi) ln(l - n,,) 



(6) 



2-q, if a + /3e^<0. 



and the g-distribution ([T]) reduces to the standard Fermi 
distribution. 

Now, the most fundamental requirement of a thermosta- 
tistical formalism is thermodynamical consistency. That 
is, the formalism must comply with the standard thermo- 
dynamical relationships among thermodynamical vari- 
ables such as entropy, energy, temperature, etc. For in- 
stance, one requires the well known relationship 



2 



between the entropy, the energy, and the temperature of 
a thermodynamical system at equihbrium to be satisfied, 
ft can be shown (see [l^, [2l| and references therein for 
details) that any thermostatistical formahsm constructed 
on the basis of the constrained externalization of an en- 
tropic functional (that is, following Jayne's maximum 
entropy prescription) complies with the thermodynam- 
ical relationships (which, in the context of Jayne's' max- 
ent formulation are usually referred to as Jaynes' rela- 
tionships). To obtain a thermodynamically consistent 
formulation one has to make the appropriate identifica- 
tions between relevant constraints and extensive thermo- 
dynamical quantities, on the one hand, and between the 
corresponding Lagrange multipliers and appropriate in- 
tensive thermodynamical quantities, on the other one. In 
the case of the formalism based upon the entropic varia- 
tional principle ([S]) the appropriate identifications are 

i 

Y.n\^, -> (8) 

i 

and 

-a//3 ^ M- (9) 

The functional (jSj, of course, is to be identified with the 
entropy of the system. It is plain from that, in order 
to compute physical quantities in a thermodynamically 
consistent way, one must not use directly the particle 
distribution given by eq. ([!]) (as is done in See 
eqs. (17-20) in [TgI) but, instead, use the effective particle 
distribution 

n\ = (l + [l + (9-l)/3(e.-Ai)]^)"', (10) 

with q defined as in ([2]). For instance, if the energy of 
an ideal Fermi gas is to comply with the basic thermody- 
namical relationship ([7|), it has to be computed (accord- 
ing to dH])) as 

£; = ^ e, (l + [1 + (g - 1) /3(e, - ^)]l^) . (11) 

i 

It must be stressed that making the identification 
does not imply any severe conceptual difficulty, since the 
fii are not probabilities and, consequently, are not nor- 
malized to unity (see (22j for a similar situation arising 
in connection with the ^-generalization of the classical 



Boltzmann distribution). In fact, it is possible to re- 
formulate the variational principle ([5]) in terms of linear 
constraints, by recourse to an appropriate re-definition 
of the entropic functional. Indeed, if one introduces the 
new entropy, 

^r = E^«("'^')' (12) 

i 

with the function Cq still defined as in Q , then the vari- 
ational principle 



. . . / 

is equivalent to the variational principle ([5]). The so- 
lution to this variational problem is given by the mean 
occupation numbers 

= (l + [l + (g-1) /?(£,- Ai)]^)"', (14) 

which are given by the same expression as the one given 
in (|10p . The thermodynamical quantities 'N and E are 
now expressed in terms of the mean occupation numbers 
(|14p using the standard linear forms. 

There is another possible modification of the maxent 
variational problem that can be implemented in order 
to recover thermodynamic consistency. One can redefine 
the entropy functional to be 

Sq^ OixY,{n,-n\)+fi^e,{n,-nX). (15) 

i i 

In this case the variation of the entropy functional with 
the standard constraints N — hi and E — e^fii is 
given by 




= (16) 



where a\ and /3i are constants as q to be determined. 
Now if the following choices, a\ = a and /3i = /3 are 
made, the solution to the variational problem given in 
eq (jl6p is formally given by eq(IT]) where a and /3 are now 
determined by the standard equations of constraint since 
both variational equations, ([5]) and (jl6p . are equivalent. 

However, in order for this to be a useful result, it is es- 
sential that thermodynamic consistency is also preserved. 

In order to verify this it is necessary to evaluate {-q^^vn ■ 
This can be done in the following manner. 
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qn. 



9-1 



qn. 



9-1 



- 1 



(17) 



as required. 

Thermodynamic consistency is therefore preserved 
when the following identifications are made 



i 

^^fiiei -J> E, 



and 



(18) 



a strictly convex function. Thus, this redefinition in- 
troduces no additional structure into the entropy density. 

In the present work we have considered the q-deformed 
quantum entropy functional given in terms of Fermi- 
Dirac distribution functions. Similar results can be ob- 
tained for quantum systems given in terms of Bose- 
Einstein distribution systems. It also is interesting to 
note that analogous results can easily be obtained for the 
q-deformed classical entropy functional. Modifying the 
entropy functional admits the possibility of using stan- 
dard linear constraints such that thermodynamic consis- 
tency is preserved. 

Summing up, our main conclusions are the following. 



/3 ^ l/fcT, 
-a/P ^. (19) 

It should also be noted that in the limit g — ?> 1, the last 
two terms in ([15]) cancel and it again reduces to ^ . Re- 
quiring S** — > in the limit T places the restriction 
q<2. 
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FIG. 1: S* IL for m = 1 MeV and <5 = lO^MeV. 

To illustrate the behavior of this redefined entropy 
functional, consider the toy model of a one-dimensional 
Fermi gas. In this case the entropy density (entropy per 
unit length) can be calculated from 




\C^{n) + (e - y)P(n - n«)].g(e)de (20) 



where (5 is a high energy cutoff introduced for numerical 
purposes, and g{e)de = -^de is the one-dimensional 
density of states per unit length. For m = IMeV 
and ,5 = 10^ MeV, FIGH] shows S*/L for various 
values of q. The entropy density increases at a given 
temperature for increasing values of q but remains 



• First, a straightforward application of the varia- 
tional principle ([5]) is not consistent, from the ther- 
modynamical point of view, with the use of the 
standard identifications N fii and E 

Cifii for the total number of particles and the 
total energy, respectively (a similar problem occurs 
with other thermodynamical quantities that in the 
standard thermostatistical formalism are expressed 
as linear functions of the mean occupation num- 
bers). 

• Thermodynamic consistency can be recovered by 
using the identifications N — J^i ^ind E 
^i^ifil- This approach does not lead to serious 
conceptual problems, because the mean occupation 
numbers are not probabilities and are not normal- 
ized to one. 

• Alternatively, thermodynamic consistency can be 
recovered by appropriately redefining the entropy 
functional. In the present work we have consid- 
ered two alternative ways of implementing this last 
procedure. 

Thermodynamical consistency is certainly a strong and 
fundamental constraint in the development of extended 
or generalized thermostatistical formalisms of physical 
significance. However, in the case of the q-deformed 
Fermi-Dirac statistics, this requirement alone does not 
determine unequivocally a unique non-extensive gener- 
alization of the standard statistics for fermions. As 
we have pointed out, there are several ways to imple- 
ment a thermostatistical formalisms for the q-deformed 
Fermi-Dirac distribution in a thermodynamically consis- 
tent way. Only when more experimental data is available, 
and more applications to concrete quantum many-body 
systems are investigated, will one definitively be able to 
ascertain which choices of the concomitant entropy func- 
tional and constraints should be used and under which 
circumstances. Any further developments along this lines 
will be very welcome. 



4 



Acknowledgments the programs No. FQA/[-2445 and No. FQM-207 of the 

Junta de Andalucia (Spain). 

We would hke to thank J. Cleymans for bringing this 
to our attention This work was partially supported by 



[1] F. Buyukkilic and D. Demirhan, Phys. Lett. A 181 

(1993) 24. 

[2] F. Buyukkilic, D. Demirhan, and A. Gulec, Phys. Lett. 

A 197 (1995) 209. 
[3] F. Pennini, A. Plastino, and A.R. Plastino, Phys. Lett A 

208 (1995) 309. 
[4] S. Curilef, Phys. Lett. A 218 (1996) 11. 
[5] S. Curilef, Zeitschrift fur Physik B 100 (1996) 433. 
[6] U. TirnakU, F. Buyukkihc, D. Demirhan, Physica A 240 

(1997) 657. 

[7] U. Tirnakh, F. Buyukkihc, D. Demirhan, Phys. Lett. A 

245 (1998) 62. 
[8] M. Ubriaco, Phys. Rev. E 60 (1999) 165. 
[9] U. Thuakh and D. Torres, Eur. Phys. J. B 14 (2000) 691. 
[10] R. Silva, D.H.A.L. Anselmo, and J.S. Alcaniz, EPL 89 
(2010) 10004. 

[11] C. Tsahis, Introduction to Nonextensive Statistical Me- 
chanics: Approaching a Complex World, Springer, New 
York, 2009. 

[12] M. Gcll-Mann and C. TsaUis (Eds.), Nonextensive En- 
tropy - Interdisciplinary Applications, Oxford University 
Press, 2004. 



[13] C. Tsallis, M. Gcll-Mann and Y. Sato, Proc. Natl. Acad. 

Sci. USA 102 (2005) 15377. 
[14] J. Rozynek and G. Wilk, Journ. Phys. G 36 (2009) 

125108. 

[15] J.M. Conroy and H.G. Miller, Phys. Rev. D 78 (2008) 
054010. 

[16] F.I.M. Pereha, R. Silva, and J.S. Alcaniz, Phys. Rev. C 

76 (2007) 015201. 
[17] H.G. Miller, F.C. Khanna, R. Tcshima, A.R. PLastino, 

and A. Plastino, Phys. Lett. A 359 (2006) 357. 
[18] A.M. Teweldeberhan, H.G. Miher, and R. Tegen, Int. J. 

Mod. Phys. E 12 (2003) 395. 
[19] H. Uys, H.G. Miller, and F.C. Khanna, Phys. Lett. A 

289 (2001) 264. 
[20] A.M. Teweldeberhan, A.R. Plastino, and H.G. Miller, 

Phys. Lett. A 343 (2005) 71. 
[21] A.R. Plastino, A. Plastino, H.G. Miller, and H. Uys, As- 
trophysics and Space Science 290 (2004) 275. 
[22] J. Cleymans, G. Hamar, P. Levai, and S. Wheaton, J. 

Phys.G: Nucl. Part. Phys. 36 (2009) 064018. 



